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00 Abstract: An impulse response matrix of a received signal in a TDMA communications system is approximated using a plu- 
^ rality of indirect variables of a linear complex vector. The indirect variables are used for synchronizing to the received signal (14A, 
^ 14B) and for tracking (16A, 16B) and frequency offset estimation (22) during successive samples of the received signal, the samples 

being equalized in dependence upon the indirect variables (18). A demodulated signal (28) is derived from the equalized received 
O signal samples. Individual synchronization (14A. 14B) and tracking units (16A, 16B), and a single equalizer (18), can be provided 

for a two-antenna receiver. Tracking errors can be used to adapt a parameter of the equalizer to reduce interference in the received 
^ signal (30). 
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CELLULAR COMMUNICATIONS SYSTEM RECEIVERS 

This invention relates to cellular coinm\inications 
system receivers, and is particularly concerned with such 
receivers for use in mobile radio systems, such as TDMA 
5 (time division multiple access) systems - 

Background 

It is well known that it is necessary in a 
receiver of a cellular mobile radio system to recover each 
communicated signal under varying and challenging 
conditions. These conditions include, for example, the 
presence of multi-path signals and fading resulting in low 
signal-to-noise ratio (SNR) , the presence of strong co- 
channel interfering (CCD signals, and Doppler effects due 
to relative movement of the signal transmitter and 
receiver, as well as typical constraints due to factors 
such as limited channel bandwidth and equipment tolerances. 
In the case of a TDMA system, it is necessary to recover 
the timing and synchronize to the time division multiplex 
(TDM) frames and time slots of a received digital 
communications signal. It is also desirable to provide the 
receiver with the least possible cost and computational 
complexity. 

International Piiblication Nimber WO 97/08867 
dated March 6, 1997, in the name of Northern Telecom 
25 Limited and entitled Timing Recovery And Frame 

Synchronization In A Cellular Communications System", 
discloses a method of timing recovery in which indirect 
variables of a linear complex vector are estimated using a 
maximum likelihood criterion in order to recover sampling 
30 delay and hence the timing and frame synchronization of the 
received signal. 
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There remains a need to provide improved 
receivers for cellular communications systems. 

Summary of the Invention 

According to one aspect of this invention there 
is provided a method of processing samples of a received 
signal to produce a demodulated signal, comprising the 
steps of representing an impulse response matrix for the 
received signal using a plurality of indirect variables of 
a linear complex vector, and synchronizing to the received 
signal samples in dependence upon the indirect variables; 
characterized by the steps of: tracking the indirect 
variables for successive received signal samples; 
equalizing the successive received signal samples in 
dependence upon the tracked indirect variables; and 
producing the demodulated signal in response to the 
equalized received signal samples. 

Thus in accordance with this method the indirect 
variables track variations in distortions, such as delay, 
fading, and phase distortions, and their use is extended to 
20 signal processing steps for producing the demodulated 

signal, thereby facilitating an improved performance of the 
entire receiver of a communications system. 

The received signal may in particular be a signal 
of a TDMA communications system, and the step of 

25 synchronizing to the received signal samples can comprise 
matched filtering the received signal samples to produce 
the plurality of indirect variables, and determining a 
maximum of a function of the indirect variables to 
determine synchronization. In a particular embodiment of 

30 the invention described below, there are four indirect 
variables and said function is a function of only two of 
the indirect variables . 
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The step of tracking the indirect variables for 
successive received signal samples can comprise recursively 
filtering initial values of the indirect variables, 
established during the synchronizing step, in dependence 
5 upon the successive received signal samples, and can also 
comprise a step of estimating frequency offset in 
dependence upon the successive received signal samples . 
This enables the tracking to be effective over time slots 
in a TDMA system operating at high frequencies, for example 
10 2-4 GHz, despite rapid changes due to Doppler effects 

arising from relative movement between a transmitter and a 
receiver of the received signal. 

The step of equalizing the successive received 
signal samples can comprise adaptively changing an 
15 equalizer parameter in dependence upon tracking errors for 
successive received signal samples to reduce co-channel 
interference with the received signal* 

The method can advantageously be applied in a 
dual antenna receiver arrangement in which said indirect 

20 variables are produced and tracked individually in respect 
of samples of a received signal from each of two spaced 
antennas, received signal samples from the two antennas 
being combined and equalized in dependence upon a 
combination of the indirect variables in respect of the two 

25 antennas . 

Another aspect of the invention provides 
apparatus for producing a demodulated signal from samples 
of a received signal, comprising a synchronization unit 
responsive to the received signal samples for producing a 
30 linear complex vector comprising a plurality of indirect 
variables having initial values corresponding to a 
synchronized state, the apparatus being characterized by: 
a tracking unit responsive to the initial values of the 
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indirect variables and to the received signal samples to 
produce tracked values of the indirect variables for 
successive received signal samples; an equalizer responsive 
to the tracked values of the indirect variables to equalize 
5 the successive received signal samples; a feedback path 
from the equalizer to the tracking unit to facilitate 
producing the tracked values of the indirect variables by 
the tracking unit; and a demodulator responsive to the 
ecjualized received signal samples to produce a demodulated 
10 signal. 

The synchronization unit can comprise a plurality 
of finite impulse response filters for matched filtering of 
the received signal samples to produce the plurality of 
indirect variables. The tracking unit can comprise a 
15 recursive filter for recursively filtering the indirect 

variables in dependence upon the successive received signal 
samples . 

The apparatus can also include a frequency offset 
estimator coupled to the tracking unit for modifying the 
20 tracking of the indirect variables in accordance with 
estimated frequency offset in dependence upon the 
successive received signal samples. 

For reducing co-channel interference, the 
apparatus advantageously includes a unit, responsive to 
25 tracking errors determined by the tracking unit for 
successive received signal samples, for estimating an 
interference correlation matrix to adaptively change a 
parameter of the equalizer. 

The apparatus can include respective 
30 synchronization and tracking units for samples of a 
received signal from each of two spaced antennas, the 
equalizer being responsive to the tracked indirect 
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variables for both antennas to combine and equalize the 
received signal samples from the two antennas. 

Brief Description of the Drawings 

The invention will be further understood from the 
5 following description by way of example with reference to 
the accompanying drawings, in which: 

Fig. 1 schematically illustrates a dual antenna 
TDMA cellular radio system receiver using indirect 
variables in accordance with an embodiment of the 
10 inventions- 
Fig. 2 schematically illustrates indirect 
variable synchronization and frame synchronization units of 
the receiver of Fig. 1; 

Fig. 3 schematically illustrates a linear 
15 indirect variable equalizer of the receiver of Fig. 1; 

Fig. 4 schematically illustrates indirect 
variable tracking, frequency offset estimation, and K<1> 
units of the receiver of Fig. 1; 

Fig, 5 schematically illustrates a hard limiter 
20 of the receiver of Fig. 1; 

Fig. 6 schematically illustrates an optional 
interference correlation matrix estimator of the receiver 
of Fig. 1; and 

Fig. 7 schematically illustrates a single antenna 
25 TDMA cellular radio system receiver using indirect 

variables in accordance with another embodiment of the 
invention. 

Detailed Description 

Embodiments of the invention are described below 
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in the context of a receiver for use in a TDMA cellular 
system compatible with EIA/TIA document IS-54-B: Cellular 
System Dual Mode Mobile Station — Base Station Compatibility 
Standard (Rev. B) and later documents referred to as IS-136 
5 and IS-136+, minimum performance standards for which are 
specified in a document referred to as IS-138. For 
brevity, this is referred to below simply as the TDMA 
system and the IS-54 standard. However, the principles of 
the invention are also applicable to other TDMA systems and 
10 to other types of communications system receiver. 

Such a TDMA system requires a receiver having a 
high performance for operation in various radio channel 
conditions which include a frame duration of 20 ms, a 
channel bandwidth of 30 kHz, a multi-path delay spread of 

15 up to 42 ps (one TDMA symbol) with equal powers of the 

multi-path signals, Doppler frequency up to 200 Hz, and the 
presence of up to 3 strong co-channel interference signals. 
The receiver is desired to provide the best possible 
sensitivity and multi-path fading reception, despite these 

20 conditions, with the least possible . computation complexity 
and cost. 

It is well known to enhance reception of radio 
channels subject to fading by using two (or more) spaced 
antennas whose respective receive path signals are combined 
25 In a desired manner, and the first embodiment of the 
invention described below relates to a dual antenna 
receiver arrangement. However, the invention is also 
applicable to a single antenna receiver, as described later 
below. 

30 Despite the use of two antennas, a design of 

receiver for operation in the various conditions outlined 
above presents a significant challenge. As indicated 
above, one such design described in p\iblication WO 97/08867 
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makes use of indirect variables to recover timing and frame 
synchronization of the received signal. 

As described in that publication, each component 
of an impulse response matrix G(t) is approximated by a 
5 linear combination, plus a constant term, of a pair of 
functions (pi(x) and 92 (t). Several examples of function 
pairs are given. This leads to introduction of a variable 
On which is a 3-dimensional complex vector constituted by 
the transpose of three indirect variables <t>i,n/ ^2,nr and 
10 (|>3,n/ which are used in the processes of timing recovery and 
frame synchronization. The recovered timing and 
synchronization parameters are then used in conventional 
manner for deriving the content of the received signals. 

The present invention also uses indirect 
15 variables, but does not merely use them for timing recovery 
and frame synchronization. Rather, the invention 
recognizes that parameters such as the sample timing and 
frame synchronization are only means to the end of 
recovering the content of the received signal, and that 
20 these parameters do not necessarily provide any value for 
themselves. Instead, the invention makes use of indirect 
variables substantially throughout the entire receiver, and 
then recovers the content of the received signal from the 
indirect variables at the demodulator. This facilitates 
25 achieving an improved receiver performance- In particular, 
as described further below, the indirect variables can be 
used in tracking channel changes, Doppler and other 
frequency offsets, providing equalization, and also in 
reducing the adverse effects of co-channel interference. 

30 Referring to Fig. 1, which illustrates a dual 

antenna TDMA system receiver using indirect variables in 
accordance with an embodiment of the invention, two spaced 
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antennas lOA and lOB are coupled to respective receiver 
front end units 12A and 12B, each of which includes a radio 
frequency receiver, dovm converter, sampler, and analog-to- 
digital converter of known form to provide at its output 

5 digital complex signal samples and i^, the subscript n 
denoting the sample number and the superscript A or B 
denoting the antenna. These signal samples are supplied to 
a respective one of two indirect variable (IN. VAR. ) 
synchronization units 14A and 14B, to a respective one of 
10 two independent variable tracking units 16A and 16B, and to 
a linear independent variable equalizer 18 of the receiver 
of Fig. 1. 

The receiver of Fig. 1 also includes, commonly 
for the signals of the two antennas, a frame 
15 synchronization unit 20, a frequency offset indirect 

variable estimator 22, a hard limiter 24, a calculation 
unit 26, and a demodulator 28. Signal connections among 
these various units of the receiver are shown in Fig. 1 and 
are further described below. In addition, the receiver can 
20 optionally include an interference correlation matrix (ICM) 
estimation unit 30 which is shown with its connections in 
dashed lines in Fig. 1. 

In order to understand the further description 
below, it is expedient to consider a mathematical 
background which leads to such understanding. This 
consideration is for a two-path signal model in an IS-54 
TDMA system, in which as is well known signals are 
communicated using n/4 -shifted DQPSK (differential 
quadrature phase shift keyed) signal symbols in non- 
overlapping time slots each of which comprises data 
symbols, known synchronization symbols (sync word), and 
known CDVCC symbols. The known symbols, in particular the 
sync word, facilitate determination and tracking of 
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parameters of the signal received during the time slot, 
these parameters including for example carrier phase which 
can vary during the time slot. 

With sampling as is usual at twice the clock 
5 frequency, a discrete observation model for the received 
signal samples has the form: 

M / ^Tjt rp\ M / ^rp rji\ 

Vi = Z "'I " - 2J - -^2 - kT - -J + T,i (1) 

where yi is the complex observation sample, i denotes the 
sample number from 1 to 2N+1 in the observed data sequence, 
10 Sk are the known complex symbols in the sync word of M 

symbols, u^ = -^^^xp^xi/q) are unknown complex amplitude- 
phase multipliers for the different independently fading 
paths m=l and m=2, the samples i on each path m having 

power ^ and the average power of each path being half the 
15 average signal power, T is the symbol spacing or clock 
frequency period, Xi and Xz are unknown delays of the two 
paths, g(t) is the impulse response of concatenated 
transmitter and receiver filters, given by: 

, X /"sin^TCt / t)Y cos(a7ct/T) ^ 

20 where a is the filter roll-off coefficient, and Tii is a 
noise sequence of complex Gaussian random variables with 
zero mean, variance 2gt^, and correlation function 2cr7^g(({i- 
j)T)/2) between two random variables Tii and T]j. 

The channel model of Equation (1) can be written 
25 in the form: 

Yi = E - - kT - + i - -^2 - - |j + T,i (2) 

where = the amplitude-phase multipliers during 

the sync word and Sj^ = S)^(s„) are transformed symbols of the 
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sync word. Assuming that the amplitude-phase multipliers 
are constant during the sync word, then Equation (2) can be 
rewritten in matrix form as: 

^ = G(Ti)SDi + g(x2)su2 + (3) 

5 where \ = [y2a-i Van - y2n+2N-2 Yan+aN-iF is a (2N+1)- 
dimensioned observation vector, 

^ = han-i T]2n 'n2n+2N-2 'n2n+2N-iF a (2N+1 ) -dimensioned 
vector of correlated noise samples, S = (sjq) [sq &l s^-i s^F 

is an (M+1) -dimensioned vector of known symbols, and G(x) 
10 is the impulse response matrix given by: 

g(- -c) g(- - t) ... g(- T - mt) 



G(x) = 



g(- T + T / 2) g(- T - T / 2) 
g(- T + t) g(- t) 

g(- T + NT) g(- x + (n - l)r) 



g(- X + T / 2 - mt) 
g(- X - (m - l>r) 

g(- X + (n - M)r) 



According to the IS-54 standard, 
Xi e (- Lpj,T / 2; Lp^T / 2) i = 1, 2 |xi - Xj] < T 
where Lpr is the number of sample spacings in an 

15 uncertainty range of delay, abbreviated below to L. With 
this uncertainty range deteinnined to be L==2 (path delays 
within two sampling intervals) , then to ensure that the 
sync word symbols are all within an observation set the 
above impulse response matrix must be increased by two 

20 initial rows and two final rows, so that with N=M it 
becomes a matrix with 2M+5 rows and M+1 columns with 
components Gij=g(-x+(i-l-L)T/2- ( j-l)T) where i is a row 
index from 1 to 2M+5 and j is a column index from 1 to M+1, 

In publication WO 97/08867 an approximation 
25 g(x + iT / 2) = ai,iq>i(x) + a2,i92W + ^3,± i = ... - 2, - 1, 0, 1, 2, ... 
or, in matrix form, g(x) = Ai9i(x) + A292W + ^3 is used for a 
single path channel to approximate the impulse response 
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matrix G(t), where ai,i;- 3L2,if ^rid 33,1 are approximation 
coefficients and Ai, Azf and A3 are approximation matrices 
with these coefficients. This approximation is reasonable 
for values of the variable x within the interval [-T/2; 
5 T/2], but is not sufficient for two paths for which, with 
frame synchronization established, the uncertainty range is 
still [-3T/2; 3T/2] . A good accuracy with this greater 
interval has been found with four terms, i.e.: 
g(T: + iT / 2) = ai,i<pi(x) + az^^q^i^) + ag^iCpgCt) + a^^^<f^{x) 

° i = ... -2, -1, 0, 1,2,... 



or g{x) s Ai<Pi(t) + A2CP2W + A393(t) + A494(t) (4) 

Various approximation functions can be used to 
provide a desired approximation accuracy, and the invention 
15 is not limited to any particular set of approximation 

functions. As one example, the approximation functions may 
be: 

(Pi(t) = sin{7rc / 2T) 92W " cos (ttc / 21) 

93(x) = sin (ttx / 4T) 94 W = cos (ttc / 4T) 

20 As another example, the approximation functions may be: 

9i('t) = g(T^) <P2('c) = g('c) 

93(x) = -1.22g(x) + g(x / 2) 94 (x) = -g(x) + 1. 41g(x / 2) 
where g(x) is the Hilbert transform of the function g{t). 
These approximation functions provide a signal-to- 

25 approximation noise ratio of about 30 dB in a range of x 
from -T to T. As can be appreciated, a preferred set of 
approximation functions, and the number of functions in the 
set, depends on the desired approximation accuracy (signal- 
to-approximation noise ratio) and directly affects the 

30 resulting complexity of implementing the approximation 
functions in the receiver. 
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Because the difference in delays between the two 
paths is not more than one symbol spacing interval. 
Equation (3) above can be rewritten, using the 
approximation functions, in the form: 

5 

If a 4-dimensioned vector On of indirect 
variables is defined by: O^ s ((t>i,n/ 4>2,nr <t>3,n^ Kn)^ 
where ^^^^ = <?^{xM + 9i('^2)Un i = 1/ 2, 3, 4 (6) 

10 are the indirect variables, then the model of Equation (5) 
can be expressed in the form: 

- BO^ + (7) 
where B = [(AiS) (AjS) (A3S) (A4S)] is a known matrix because Ai to 

A4 comprise fixed coefficients and S is the known sync 
15 word. 

In this context, in the receiver of Fig. 1 the 
indirect variable synchronization units 14A and 14B serve 

to produce initial values Oq and Og respectively of the 

indirect variable vector On for synchronization, and the 
20 indirect variable tracking units 16A and 16B serve, in 
conjunction with the frequency offset indirect variable 
estimation unit 22, to track the indirect variable vector 

On throughout a time slot to produce tracked values O^ and 

On respectively. The linear indirect variable equalizer 
25 18 comprises a Kalman filter which is controlled by the 
tracked values O^ and O^ of the indirect variable vector 
to combine and recursively filter the received signal 
samples y^ and y^, thereby producing a received and 

SUBSTITUTE SHEET (RULE 26) 



wo 01/89165 PCTmUOO/00181 

13 

equalized signal vector Se,n- This vector is limited by 
the hard limiter 24, from the output of which the KO 
calculating unit produces a feedback control signal K<i>,n 
for the tracking units and the demodulator 28 produces a 
5 demodulated signal on an output line 32 of the receiver. 
The various units of the receiver are further described 
below . 

Fig. 2 illustrates one form of the indirect 
variable synchronization unit 14A (the unit 14B is similar) 
10 and one form of the frame synchronization unit 20, for 

producing the initial value Oq for synchronization. The 
unit 14A comprises four finite impulse response (FIR) 
filters (FIR-1 to FIR-4) 40 which are supplied with the 

received signal samples , a calculation unit 42, a delay 
15 line 44 comprising delay elements each providing a delay of 
one symbol spacing interval T, and a selector 46. The 
frame synchronization unit 20 comprises a combiner 48, a 
delay line 50, and a maximum detector 52. 

In order to simplify matrix inversion, the model 
20 of Equation (7) is divided into even and odd sample sets so 
that the model can be 'expressed as: 

and, because B is a known matrix, the indirect variable 
vector can be initially determined by a matched filtering 
25 represented by: 

^odd _ (Bodd^odd)" ^odd^^ 



,even 



~ V^even^even/ "even-4i 
in which the first three terms on the right-hand side of 
each equation can be pre-calculated and stored. It is 
30 observed that here and below the various signal processing 
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operations produce results which are estimates rather than 
the precise values of the respective signals. In Fig. 2, 
the FIR filters 40 perform this matched filtering fvmction^ 
the filters FIR-1 and FIR-3 corresponding to the first and 

5 third rows respectively of the matrix (B^dSo^id)"^Sodd^^*^ 

the filters FIR-2 and FIR-4 corresponding to the second and 

fourth rows respectively of the matrix (Beven^even^^BevenYn 



^even 



Consequently, the outputs of the FIR filters 40, 
which are supplied to inputs of the delay line 44, 

10 constitute the indirect variable vector On in accordance 
with the above model, but its synchronization, i.e. the 
value n which provides a reference timing point, is not yet 
determined. This is deteimined as described below using, 
for simplicity, only the first two of the approximation 

15 functions described above, the outputs of only the filters 
FIR-1 and FIR-2 being supplied to the calculation unit 42, 
which calculates and produces at its output a value modOn. 
The calculation carried out by the unit 42 is dependent 
upon the particular approximation functions which are used 

20 as described above. For the functions <Pi(x) = g(x) and 
92 W = ^K'^) referred to^ above, for example. 



mod 



2^ 



where = — real] a tan 



and the calculation unit 42 



calculates modOn accordingly. 

25 The values modO^ and thus produced for the 

two paths A and B are combined by the signal combiner 48 in 
the frame synchronization unit 20 to provide for frame 
alignment of the signals from the two antennas, and the 
output of the signal combiner is supplied to the delay line 
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50 having a length corresponding to the observation window 
of the received signal samples. The maximiiin detector 52 
determines the va2jae n corresponding to a maximxim one of 
the outputs of the delay line 50, thereby determining 
5 synchronization, and supplies this to the selector 46 in 
the unit 14A. The selector 46 selects from the delay line 
44 the corresponding indirect variable vector On and 
supplies this to its output as the initial indirect 
variable vector Oq - 

10 One form of the linear indirect variable 

equalizer 18 is illustrated in Fig. 3 and comprises a 
Kalman filter including signal combiners 60 and 62, 
multipliers 64 and 66, and a delay element 68, and an 
arrangement for determining filter parameters including 

15 calculation units 70, 72, 74, and 76 and a delay element 
78. The arrangement and operation of the equalizer will be 
further understood from the following description. 

With prediction of estimates to an n-th step at 
the step n-L, the Kalman filter in Fig. 3 can be seen to 
20 provide its output Se,n in accordance with: 

S0,n = ^n^n-l^ + Kg^nG^n " / n-LXAnSn-l)) (8) 

where Ks,n is a Kalman filter gain given by: 

Ks,n = / h-iGK / n^L)l(G(«>^ / n^lJh^n / n-lG^ / h-l)] + V^^)"^ 

Vs,n/n-l = A^Vg^n-lA^ + V^,^ (9) 

where the terms of these equations can be understood from 
the following description. 

Extending a single antenna representation for the 
case of two antennas and hence two sampled signals. 
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Yn - 




Sn + 



Tin 



where the last matrix represents equivalent noise including 
indirect variable errors. This noise has the covariance 



matrix: 



V = 2^?^ -^^n^^i* 02x2 



^23c2 



5 R. 



1 g(T/2)1 

[g(T/2) 1 J 



where ^^^n/n-h ~ ^o^n-L + Qa>^ is the covariance matrix of 

the indirect variable prediction error for L+1 steps, Qo is 
the covariance matrix of exciting noise of the indirect 



10 variable vector, b(s) = 



r^i . X' L and A2 



are 



4 x (2m+l) -dimensioned approximation coefficient matrices 
as described above, and m is a number of adjacent symbols 
taken into account. A vector condition for TDMA symbols 
can be described by the equation = A^Sn-i + §n where An 
15 is a shift matrix (or a transition matrix during the CDVCC) 
and ^ is noise with covariance matrix V^^j^ = 2Q^,n which is 

a zero matrix during the CDVCC, thus when n is not in the 
CDVCC: 
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Thus referring again to Fig. 3, the calculation 
unit 70 produces G(On/n-L)/ used in the first and third of 
Equations (9) and in Equation (8), from the indirect 

variable vectors and <I>^ tracked as described below, 
5 and the calculation units 76, 74, and 72 produce the values 
of respectively the first, second, and third of the 
Equations (9) . V^^n supplied to the unit 74 is 
predetermined as indicated above, and VT|,n supplied to the 
unit 76 can also be fixed and predetermined or, as later 
10 described below, can be adaptively changed. The Kalman 
filter gain Ks^n is used by the multiplier 66 to produce, 
with the other elements 60, 62, 64, and 68 of the Kalman 
filter, the equalizer output vector Se,n from the signal 
samples y^ and y^ in accordance with Equation (8) . 

15 Referring now to Fig. 4, one form of the indirect 

variable tracking unit 16A is illustrated, the unit 16B 
being similar and connections to it also being indicated in 
Fig. 4, together with one form of the frequency offset 
indirect variable estimation unit 22 and of the KO 

20 calculation unit 26. To some extent these units make use 
of similar calculations so that they are closely inter- 
related and they are accordingly described together below. 
As illustrated in Fig. 4, the indirect variable tracking 
unit 16A comprises signal combiners 80 and 82, multipliers 

25 84, 86, 88, and 90, and a delay element 92. Although not 
shown separately in Fig. 4, the unit 16A is also supplied 

with the initial value Oq of the indirect variable vector 

On for synchronization by the unit 14A as described above. 
The frequency offset indirect variable estimation unit 22 
30 comprises signal combiners 94, 96, and 98, multipliers 100 
and 102, a delay element 104, and calculation units 106, 
108, and 110. The calculation unit 26 comprises signal 



SUBSTITUTE SHEET (RULE 26) 



wo 01/89165 PCT/RUOO/00181 

18 

combiners 112 and 114/ a delay element 116, and calculation 
units 118, 120, 122, and 124. 

The effect of frequency offset in conjunction 
with indirect variable error ^ can be expressed by: 
5 . dF,Oti + ^tn 

= dFX-l + ^f,n (10) 

where dF^ « exp(i27i^fT) and fot is frequency offset. With an 
observation model of the form; 
10 = b(s>^ + t,* = B(SJ^I + Ti« 

then if dFn-i is estimated very accurately, substituting the 
estimates for the actual variable in Equation (10) produces 
a model which does not depend upon this variable: 

15 = Ol^i +^^,n (10) 

and results in a filtering algorithm: 

= dF^_i0^.i + K^^nif/n " B(SnK-xdF^.J (11) 

where K^^^ = ^<j,,n/ n-^A^J%^n 

20 %X = ksnKn/n-lBCSn/ + 2a^R^Y (12) 

It can be appreciated here that the same Kalman gain matrix 
K<D,ii is used for signals from the two antennas, thereby 
25 simplifying the filtering algorithm and making its 
complexity independent of the number of antennas. 

It can be seen from Fig. 4 that the elements 80 
to 92 of the tracking unit 16A are arranged to implement 
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Equation (11) for the signal samples y^, the parameter dFn- 
1 being supplied from the output of the delay element 104 
and the parameter K^^n being supplied from the output of 
the unit 124, to produce the tracked indirect variable 

5 vector at the output of the signal combiner 82, The 
four Equations (12) are implemented in Fig, 4 respectively 
by the calculation units 124, 122, and 120 and the elements 
112 to 116, 

For frequency offset estimation, the following 
10 approximated model is derived from the above: 

dF^ = dF^.i + ^dF,n 



y^ = B(s^)0^.idF,_i + V 



where and are equivalent observation noise which 
15 take into account estimation errors of the indirect 
variable vectors 0^_i and <t>n-i and have the same 
covariance matrix defined by: 

Then an indirect variable filtering algorithm for 
20 dFn can be written in 'the form: 

Where K^F,a = VdE.,n<I^n-lB(Sn) V"], 

K^F.n = VdF,^0^.,'B(sJv-^^ (14) 

25 ^dF,n/n-l ^ ^dF,n-l + Q^dF 

Although Equations (14) appear to be very 
complex, they can be implemented with low complexity 
because a large number of steps are already otherwise 
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performed. For example, the inverse matrix V""^ is used 

for calculation of Ko,nr the differences in Equation (13) 
are the same as those in Equation (11) r and multiplications 

such as B(Sn)^>n_i have been used for difference 
5 calculations. 

It can be seen from Fig. 4 that the four 
Equations (14) are implemented by the calculation units 
108, 110, and 106 and the signal combiner 94 respectively. 
The multipliers 100 and 102, signal combiners 96 and 98, 
10 and delay element 104 implement Equation (13) • 

Referring to Fig. 5, the hard limiter 24 is 
supplied with the vector Se,n from the output of the 
equalizer 18, and supplies the elements Sn+m to Sn-m of this 
vector via respective ones of 2m+l stages 130 to produce 

15 hard limited elements which, with delayed versions thereof 
produced by delay elements 132, constitute the elements of 
the vector Sn that is supplied to the KO calculation unit 
26, and in particular to the unit 118 as shown in Fig. 4. 
In addition, the hard limiter 24 provides an output for the 

20 hard limited version of the element Sn-m to the demodulator 
28, which operates in- a well-known manner for demodulating 
the 7c/4-shifted DQPSK signal. Each of the stages 130, as 
represented in Fig. 5 for two such stages, provides signal 
phase rotation, hard limiting, and derotation following the 

25 7c/4-shifted DQPSK modulation rules. 

In order to optimize parameters for operation of 
the receiver as described above, particular values can be 
selected. For example, the number 2m+l of simultaneously 
estimated symbols in the equalizer 18 can be selected as 
\ 30 being 5 with m=2, and with L=2 as already indicated the 

niamber of symbols m+l+L used in the tracking units 16A and 
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16B is also 5. The integration interval for the 
synchronization units can be and in Fig. 4 Q^,o can be 
(diag(10"^*[2 1 0.5 0.5]), and Q^,dF can be 2.5*10"^ In 
addition, a fixed value of SNR, for example 17 dB,. can be 
5 used for synthesizing the above algorithms, as the actual 
sigrial-to-noise and interference ratio may be unknown. 

The receiver as described above is intended to 
provide an optimum performance in the presence of noise . 
However, in the presence of co-channel interference, the 

10 performance of the receiver can be degraded. In order to 
reduce or avoid such degradation, the receiver can also 
include the interference correlation matrix (ICM) 
estimation unit 30 shown in dashed lines in Fig. 1. This 
provides an adaptive control of the matrix VT,,n which is 

15 supplied to the linear indirect variable equalizer 18, so 
that interference cancellation is also achieved by the 
operation of the equalizer as described above. 

An analysis can be carried out in a similar 
manner to that described above in relation to the operation 
20 of the tracking units , but in respect of the information 
symbols Si in the TDMA time slot rather than the 
synchronization symbols, from which it can be determined 
that the noise covariance matrix Vi^,n provides interference 
cancellation based on differences of estimated values 

25 = yf - G{of and = yf - (^i^f^x which are already 

determined (Equation (8) above, using a slightly different 
notation) in the operation of the Kalman filter as 
described above. As illustrated in Fig. 1, these 
differences are supplied from the tracking units 16A and 
30 16B to the ICM estimation unit 30 to enable estimation of 
\ the matrix V|,,n- This matrix is determined by: 
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= 2R^ = 21 



where a| and <t| are unknown variances (real variables), rAs 
is an unknown correlation coefficient (complex variable), 
and R-q is the known covariance matrix already specified 
5 above - 

Fig. 6 illustrates one form of the ICM estimation 
unit 30, in which the unknowns a^, a^, and r^s are each 
averaged over a desired number NA of samples by respective 
delay lines 140 and summing units 142, the outputs of which 

10 are supplied to a calculation unit 144 to determine the 
covariance matrix V,| in accordance with the above equation. 
The unit 30 also includes multipliers 146, 148, 150, 152, 
and 154, and transpose units 156 and 158, which serve to 
produce the elements of the matrix as inputs to the delay 

15 lines 140 for averaging. Thus it can be seen that this 
adaptive operation of the receiver adds very little 
complexity to the receiver, but can substantially improve 
the performance of the receiver in the presence of co- 
channel interference. 

20 Although a dual antenna embodiment of the 

invention has been described above in detail, it should be 
appreciated that the invention is not limited in this 
respect, and it may also be applied to a single antenna 
receiver as illustrated in Fig. 7. Thus as shown in Fig. 

25 7, a single antenna 160 is coupled via a receiver front end 
unit 162 whose output signal samples yn are supplied to an 
indirect variable synchronization unit 164, an indirect 
variable tracking unit 166, and a linear indirect variable 
equalizer 168. The synchronization unit 164 provides an 

30 initial synchronization vector Oo to the tracking unit 166, 
no frame alignment being required because in this case 
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there is only one signal path. The tracking unit 166 
provides a tracked vector O^r and this is equalized by the 
equalizer 168 to produce a resulting signal for 
demodulation by a demodulator 170. It can be appreciated 
5 that the units 164, 166, and 168 in the receiver of Fig. 1 
can -use similar principles to those described above in 
order to provide improved single-antenna receiver 
performance through the use of indirect variables for all 
of the signal processing in the receiver prior to the 
1 0 demodulator . 

It can be appreciated that although as described 
above the frequency offset estimation unit 22 is provided 
as is preferred to compensate for frequency of f sets, which 
may be due to local oscillator frequency variations and, 
15 especially, due to Doppler effects, in other embodiments of 
the invention this unit can be omitted. 

In addition, although the description above 
refers to, and the drawings illustrate, particular units 
such as calculation units, signal combiners, multipliers, 
20 delay elements, and so on, it should be appreciated that in 
practice the functions of all of these units can 
conveniently be carried out by one or more digital signal 
processors or application-specific integrated circuits. 

Thus although particular embodiments of the 
25 invention have been described above, it can be appreciated 
that these and numerous other modifications, variations, 
and adaptations may be made without departing from the 
scope of the invention as defined in the claims. 



SUBSTITUTE SHEET (RULE 26) 



wo 01/89165 PCT/RUOO/00181 

24 

CLAIMS: 

1. A method of processing samples of a received 

signal to produce a demodulated signal, comprising the 
steps of representing an impulse response matrix for the 
5 received signal using a plurality of indirect variables of 
a linear complex vector, and synchronizing to the received 
signal samples in dependence upon the indirect variables; 

characterized by the steps of: 

tracking the indirect variables for successive 
10 received signal samples; 

equalizing the successive received signal samples 
in dependence upon the tracked indirect variables; and 

producing the demodulated signal in response to 
the equalized received signal samples. 

15 2. A method as claimed in claim 1 wherein the 

received signal is a signal of a TDMA communications 
system. 

3. A method as claimed in claim 2 wherein the step 
of synchronizing to the received signal samples comprises 

20 matched filtering the received signal samples to produce 
the plurality of indirect variables, and determining a _ 
maximum of a function of the indirect variables to 
determine synchronization. 

4. A method as claimed in claim 3 wherein there are 
25 four indirect variables and said function is a function of 

only two of the indirect variables - 

5. A method as claimed in any of claims 1 to 4 
wherein the step of tracking the indirect variables for 
successive received signal samples comprises recursively 

30 filtering initial values of the indirect variables. 
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established during the synchronizing step, in dependence 
upon the successive received signal samples, 

6. A method as claimed in any of claims 1 to 5 
wherein the step of tracking the indirect variables for 

5 successive received signal samples comprises a step of 
estimating frequency offset in dependence upon the 
successive received signal samples. 

7. A method as claimed in any of claims 1 to 6 

• wherein the step of equalizing the successive received 
10 signal samples comprises adaptively changing an equalizer 
parameter in dependence upon tracking errors for successive 
received signal samples to reduce co-channel interference 
in the received signal. 

8 . A method as claimed in any of claims 1 to 7 

15 wherein said indirect variables are produced and tracked 
individually in respect of samples of a received signal 
from each of two spaced antennas, and received signal 
samples from the two antennas are combined and equalized in 
dependence upon a combination of the indirect variables in 

20 respect of the two antennas. 

9. Apparatus for producing a demodulated signal from 
samples of a received signal, comprising a synchronization 
unit responsive to the received signal samples for 
producing a linear complex vector comprising a plurality of 

25 indirect variables having initial values corresponding to a 
synchronized state, the apparatus being characterized by: 

a tracking unit responsive to the initial values 
of the indirect variables and to the received signal 
samples to produce tracked values of the indirect variables 
30 for successive received signal samples; 

an equalizer responsive to the tracked values of 
the indirect variables to equalize successive received 
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signal samples; 

a feedback path from the equalizer to the 
tracking umit to facilitate producing the tracked values of 
the indirect variables by the tracking unit; and 
5 a demodulator responsive to the equalized 

received signal samples to produce a demodulated signal. 

10. Apparatus as claimed in claim 9 wherein the 
synchronization unit comprises a plurality of finite 
impulse response filters for matched filtering of the 

10 received signal samples to produce the plurality of 
indirect variables. 

11. Apparatus, as claimed in claim 9 or 10 wherein the 
tracking unit comprises a recursive filter for recursively 
filtering the indirect variables in dependence upon the 

15 successive received signal samples. 

12. Apparatus as claimed in any of claims 9 to 11 and 
including a frequency" of f set estimator coupled to the 
tracking unit for modifying the tracking of the indirect 
variables in accordance with estimated frequency offset in 

20 dependence upon the successive received signal samples. 

13. Apparatus as claimed in any of claims 9 to 12 and 
including a unit, responsive to tracking errors determined 
by the tracking unit for successive received signal 
sait^leS/ for estimating an interference correlation matrix 

25 to adaptively change a parameter of the equalizer to reduce 
co-channel interference in the received signal. 

14. Apparatus as claimed in any of claims 9 to 13 and 
including respective synchronization and tracking units for 
samples of a received signal from each of two spaced 

30 antennas f wherein the equalizer is responsive to the 
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tracked indirect variables for both antennas to combine 
and equalize the received signal samples from the two 
antennas . 
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